On the lifetime of the 2\ state in 10 C 
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The lifetime of the J' r =2+ state in 10 C was measured using the Doppler Shift Attenuation Method 
following the inverse kinematics p( 10 B,n) 10 C reaction at 95 MeV. The 2+ state, at 3354 keV, has 



t = 219±(7) stat ±(10) s „ a fs corresponding to a B(E2) | of 8.8(3) e 2 



This measurement, 



combined with that recently determined for 10 Be (9.2(3) e 2 fm 4 ), provides a unique challenge to ab- 
initio calculations, testing the structure of these states, including the isospin symmetry of the wave 
functions. Quantum Monte Carlo calculations using realistic two- and three-nucleon Hamiltonians 
that reproduce the 10 Be B(E2) value generally predict a larger 10 C B(E2) probability but with 
considerable sensitivity to the admixture of different spatial symmetry components in the wave 
functions, and to the three-nucleon potential used. 
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I. INTRODUCTION 

A new generation of ab-initio calculations based on 
realistic nucleon-nucleon forces has deepened our under- 
standing of how nuclei work. Comparison of their predic- 
tions to precise new experimental data [H-Q has guided 
improvements in the calculations, both in the computa- 
tional methods used and in the underlying Hamiltonians. 
Electromagnetic decay rates have already proven to be 
surprisingly sensitive for probing 3-body forces Q. In 
this paper we examine mirror symmetry in nuclei, com- 
paring 10 Be to 10 C, the lightest T = 1 mirror pair for 
which bound excited states exist. We have precisely mea- 
sured the B(E2;2^ — > Of) from the only bound state in 
10 C and compared it to its analog in 10 Be to probe the 
subtleties of the symmetry. The experimental result is 
difficult to reproduce using our Green's Function Monte 
Carlo (GFMC) method. A variety of calculations have 
been performed in order to better understand the issues. 

Conceptually, the A = 10 mirror nuclei are interest- 
ing. 4°Be6 can be thought of as two alpha particles with 
isospin T = 0, bound to a correlated pair of s-wave neu- 
trons, which are usually outside the alpha clusters. Simi- 
larly 6 C4 consists of two alpha particles surrounded by a 
loosely bound pair of protons. Naively one might expect 
the two protons to result in a bigger B(E2) value for 10 C 
than for 10 Be. For example, a simple classical isoscalar 
liquid drop model @ is driven by charge: the carbon 
isotope has a larger quadrupole moment (~ Z) and the 
decay strength is expected to be larger in 10 C compared 
to 10 Be by a factor (6/4) 2 = 2.25. In contrast, sim- 
ple shell models @, 0] always have properly constructed 
quantum-mechanical wave functions for states and allow 
the electromagnetic decay to be separated into isoscalar 
and isovector components. For the A = 10 system, the 
isovector contributions are predicted to be small, sym- 
metric, and act to enhance the decay in 10 Be and sup- 
press it in 10 C, resulting in a 10 C decay that should be 



~10% lower than that in 10 Be. The relative B(E2) val- 
ues in 10 C and 10 Be thus represent an interesting test of 
nuclear modeling and of the isospin dependence of elec- 
tromagnetic decays. 



Within the A = 10 system, 10 C is the more exotic part- 
ner. It has only one bound excited state, with J 71 ' = 2 + , 
at 3354 keV. It becomes unbound at 4006 keV at which 
point it can disintegrate into 9 B+p. Its mirror partner, 
10 Be, has six bound states below a breakup threshold of 
6812 keV where the 9 Be+n channel opens. The first ex- 
cited state of 10 Be also has J* = 2 + and lies at 3368 
keV, a sign that, in excitation energy at least, these con- 
figurations arc similar, despite the difference in binding 
energy. 



A pioneering Doppler Shift Attenuation Method 
(DSAM) experiment by Fisher et al, |[ in 1968 was 
aimed at understanding these issues and testing the 
intermediate-coupling shell model predictions by measur- 
ing the decay rates from the first excited states in 10 C and 
10 Be. However, the DSAM technique was new and only 
^20% precision could be achieved. 10 C was found to be 
slightly more collective, but not in glaring disagreement 
with the shell model, given the large experimental uncer- 
tainties. We have recently remeasured the lifetime of the 
2^ state in 10 Be H and determined B(E2; 2f -)■ Of ) 
9.2(3) e 2 fm 4 . Comparing this with the value B(E2; 
— > Of) = 12.2±1.9 e 2 fm 4 measured by Fisher et al, 
for 10 C still supports a larger B(E2) rate in 10 C. How- 
ever, the carbon value has substantial experimental un- 
certainties, so now, with far superior experimental tools 
and much refined theory, we can re-address this interest- 
ing problem at a level of precision which should provide 
stringent tests of ab-initio calculations. 
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II. EXPERIMENT 

The 2^ state in 10 C was populated in the inverse kine- 
matics p( 10 B, n) 10 C reaction. Beams of 10 B ions of ~ 
lpnA and 95 MeV were produced by the ATLAS ac- 
celerator at Argonne National Laboratory. Targets con- 
sisted of thin layers of CH2 on thick backings of copper 
and gold. 10 C nuclei recoiling along the beam direction 
were selected by the Argonne Fragment Mass Analyzer 
(FMA) [1| positioned 90 cm downstream of the target 
and subtending 1° around the beam direction. The 10 C 
nuclei were produced at recoil velocities of /3 = v/c ~ 
13% and emerged from the backing target layer with 
(3 ~ 10%. To satisfy the FMA energy acceptance win- 
dow, the recoils had to be further slowed down to /3 ~ 
8%. This was achieved through a series of titanium de- 
grader foils placed at the entrance to the FMA. 10 C ions 
with q = 6 + were transported to the focal plane while 
most non-interacting beam particles were rejected by the 
FMA. The selection of q = 6 + ions was very effective 
for suppressing scattered beam particles. The transmit- 
ted ions first passed through two PPAC detectors before 
being stopped 50 cm behind the focal plane in a 30-cm 
deep, two-electrode ionization-chamber operated at 50 
torr. Gamma rays were detected with the Gammasphcrc 
array [lj| consisting of 100 Compton suppressed HPGe 
detectors in 16 azimuthally symmetric rings from 9 = 34° 
to 163° relative to the beam direction. 

Figure [IJa) shows a typical energy loss (AE) versus 
total energy (Etot) spectrum obtained from the ioniza- 
tion chamber. The locus with the largest AE (solid red 
circled region) corresponds to direct population of the 2f 
state in 10 C. The large spread in total energy for these 
recoils stems from energy scattering and straggling in the 
backing and degrader foils. The wide strip of counts be- 
low the 10 C recoils is identified with 10 B scattered beam. 
The pressure in the ionization chamber was not sufficient 
to fully stop the highest energy 10 C recoils, resulting in a 
wrap around feature (punch through) in the AE versus 
E tot plot. Direct population of the ground state of 10 C 
is also observed, although due to the punch through, it 
falls appears in the same AE versus E tot area as the 10 B 
scattered beam. The 7-ray spectrum obtained by gating 
on the excited 10 C recoils is given in Fig. HJb), showing 
only the 3354-keV, 2+ — » 0+ transition in 10 C. 

The current setup offers several advantages over the 
prior DSAM measurement performed with a regular- 
kinematics reaction. To compare the present technique 
with a normal kinematics DSAM measurement, the orig- 
inal Fisher et al., experiment was also repeated with 
Gammasphere. This regular-kinematics experiment was 
performed with a 9.5-McV proton beam incident on a 
10 B target followed by a gold backing. Gamma rays were 
detected with the Gammasphere array. A comparison 
of the spectra obtained from the regular- and inverse- 
kinematics reactions is given in Fig. [5] for the same angle 
group (6 = 130°) in Gammasphere. The regular kinemat- 
ics spectrum (Fig.^a)) is very complicated, with consid- 
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FIG. 1: (Color online) (a) Energy loss versus total energy plot 
of data from the ionization chamber behind the focal plane of 
the FMA. The solid red circled region is direct population of 
the 2j" state in 10 C. Data were obtained with a 116-/ig/cm 2 
CH2 target on a 23.8-mg/cm 2 copper backing, (b) Gamma- 
ray spectrum obtained by gating on the excited 10 C recoils. 



erable background from reactions of the high-energy pro- 
tons on the 10 B target, the gold backing and scattering 
in the target chamber. In contrast, requiring detection of 
a recoiling 10 C residue by the FMA almost entirely sup- 
presses the background (Fig.[2jb)), providing a spectrum 
where only the 2^ — > 0^ transition in 10 C is observed. 
Selection of the 10 C recoils with the FMA also provides a 
well-defined angle between the recoil velocity vector and 
the direction (and subsequent detection) of the 7-rays. 
Recoil detection comes at the cost of poorer statistics; 
however, with almost no background the peak centroids 
can be reliably determined to 1-2 keV. With the regular- 
kinematics reaction, the recoiling 10 C ions have a kinetic 
energy of a few MeV, an energy regime where stopping 
powers are poorly determined and the 7-ray energy shifts 
arc only a few keV. In the inverse-kinematics reaction, 
10 C nuclei were produced at a very high recoil energy 
(E ~ 80 MeV). This allows the measurement to be per- 
formed in a velocity regime where the stopping is 99.99% 
electronic and most precisely known. The high recoil ve- 
locity also produces large Doppler shifts. With the set 
of Gammasphere angles, the range of forward-shifted to 
backward-shifted peaks spanned more than 800 keV. 
In a DSAM measurement the lifetime is derived from 
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FIG. 2: Comparison of 7-ray spectra for 10 C decays from 
DSAM in (a) regular and (b) inverse kinematics. Both spectra 
correspond to the = 130° angle group in Gammasphere. 
Reaction details are included in the figure. 



the difference between production and emission veloci- 
ties. For the inverse-kinematics reaction, the distribution 
of recoil products is very forward peaked and there are 
two solutions for residues moving along the initial beam 
direction, depending on the direction of the emitted neu- 
tron. For this particular experiment, recoils emitted at 0° 
relative to the beam direction can have recoil velocities of 
75 or 82 MeV. The FMA was always set to transmit the 
higher-energy recoil group. The initial (3 value at produc- 
tion is first measured using a self-supporting CH 2 target 
and correcting for the small energy loss in the target. 
Degrader foils of the same thickness as the backed-target 
experiments were used to replicate the DSAM measure- 
ments and ensure that the FMA entrance conditions were 
the same for both the self-supporting and backed targets. 

To determine the average velocity of the recoils at the 
time of 7-ray emission, the centroid of the 3354-kcV, 
2^ — > Of transition was determined for each of the 16 
Gammasphere angle rings. In Fig. EJa), the measured 
centroid is plotted as a function of cos(#) for a 105- 
fig/cm 2 CH2 target on a 23-mg/cm 2 copper backing. A 
fit to these data yields an average f3 at the time of 7-ray 
emission of /? — 0.12422(24). For reference, a lifetime of t 
~ would give /3 ~ 0.131 (the production velocity) while 
an infinitely long lifetime would yield /3 ~ 0.099 (velocity 
after emerging from the backing layer). The measured 
centroids are compared to the rclativistic Dopplcr shift 
formula using the best fit j3 in Fig. [2ta). Figure EJb) 
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Best fit : p = 0.12422(24) 
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FIG. 3: (Color online) (a) Measured centroids of the 3354- 
keV transition for each Gammasphere angle group. The solid 
curve shows the result of the relativistic Doppler shift formula 
taking the best-fit value of /?. (b) Similar to (a) but normal- 
ized to yl — ft 2 /(l-/3cos(0)). Lines include the best fit value 
of f3 (solid line), maximum /3 value allowed by the reaction 
kinematics (dashed line), and the @ value corresponding to 
the previous lifetime measurement (dotted line). 



illustrates the quality of the fit more clearly, by divid- 
ing the measured centroids by the function *J\ — /3 2 /(T 
/3cos(0)). Included are lines for j3 values for the best fit ft 
(f3 = 0.12422, corresponding to r = 224 fs), production (/? 
= 0.131, corresponding to r = 0), and that which would 
correspond to the previous lifetime value (f3 = 0.1265 for 
r = 154 fs). 

To determine the lifetime of the level of interest from 
the measured mean (3 value, the thicknesses of the target 
and backing layers must be known, as this defines the 
relevant transit time scale. The backing foils were pre- 
pared by rolling gold and copper foils to the appropriate 
thickness. Targets were then prepared by dipping the 
backing foils into solutions of xylene and C2H2, and sub- 
sequently evaporating the xylene. The thin layer of CH2 
on a very thick backing of gold or copper made a precise 
measurement of the CH2 thickness by a traditional a- 
gauging method impossible. The CH2 thickness was de- 
termined by comparing the yield (number of 10 C recoils 
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TABLE I: Mean lifetimes from different target and backing 
combinations determined for the 3354-keV level in 10 C. The 
CH2 thicknesses are the average of the measured values at the 
beginning and end of each DSAM measurement. 



Target (CH 2 ) 
Og/cm 2 ) 


Backing 




ATstat 


(mg/cm 2 ) 


(fe) 


(fs) 


105 


23.0 Cu 


224 


±8 


150 


31.0 Au 


215 


±10 


170 


23.8 Cu 


219 


±13 


80 


24.0 Au 


198 


±18 


300 


14.9 Au 


216 


±23 



detected in the ionization chamber, per beam intensity, 
per time) of a backed target relative to a commercially 
made self-supporting CH2 target of known thickness. A 
series of degrader foils, identical to the target backing, 
were placed behind the self-supporting CH 2 target to 
achieve the same energy recoils into the FMA as with 
the backed targets. The thickness of the CH2 layer on 
the backed targets was measured before and after the 
DSAM measurement via a relative yield measurement as 
described above. An approximately 20% reduction in 
yield was observed. For the DSAM analysis, the aver- 
age of the thickness before and after the measurement 
was used and the full range of thickness taken as the un- 
certainty As the target layers are thin, do not provide 
much stopping and the ions move very quickly through 
them, the target thickness does not contribute greatly to 
the systematic uncertainty. The characterization of car- 
bon and boron ions slo wing in CH2 , copper and gold was 



taken from the SRIM and MSTAR |12j packages. 
The two models differ in stopping powers on the order 
of 3% in the relevant velocity regime. These differences 
were incorporated into the systematic uncertainty. 

The lifetime of the 3354-keV level was measured in 
five separate experiments. The target characteristics and 
the extracted lifetimes are summarized in Table I. The 
weighted mean value is r= 219±(5) s tat ±(10) S y S tfs which 
implies B(E2; 2+ -» 0+) = 8.8(3) e 2 fm 4 . This lifetime is 
substantially longer than the previous value obtained by 
Fisher et ai, Figure U gives a comparison of the cur- 
rent measurements (solid symbols) and those of Ref. [8[ 
(open symbols). Clearly, the data are now much better 
constrained for investigating the symmetry of the wave 
functions, but theoretical guidance is needed to infer the 
meaning of the result. 



III. THEORY 

Empirically, assuming charge symmetry for the wave 
functions, the transition strengths can be written as 



B(E2) = [M(E2)} 2 /5 = [AT ■ 



BT 7 ? 



(1) 



where M(E2) is the reduced matrix element and we use 
the convention that T z = ±i for the neutron. The new 
10 Be and 10 C data can be used to infer that the isoscalar 



term is dominant, A = 3.00(1) efm 2 , while the isovec- 
tor term is much smaller, B = 0.03(3) efm 2 , a 1% ef- 
fect. In conventional shell-model calculations, isospin en- 
hancements or effective charges, e(T) are introduced to 
account for effects such as core polarization: A = ^4'e(0) 
and B = B'e(l), where A' and B' are constants derived 
for a particular wave function. Very early shell model 
calculations of Cohen and Kurath d [H for the A = 10 
system gave predictions for B(E2) strengths in the form 
ofEq. (1). These are included in Fig. S] (solid black line). 
These p-shell, mirror- symmetric wave functions provide 
the correct slope for describing the transition strengths 
between 10 C and 10 Be, however, overestimate the overall 
magnitude (Fig. [J} due to the use of very simple isoscalar 
and isovector enhancements, e(0) = 2 and e(l) = 1. Us- 
ing isoscalar and isovector enhancements now broadly ac- 
cepted for p-shell calculations [I3.[l5| . e(0) = 1.7 and e(l) 
= 0.6, one obtains 10 Bc B(E2; 2+ -> 0+) = 9.7 e 2 fm 4 
and 10 C B(E2; 2+ ->■ 0+ ) = 9.1 e 2 fm 4 , close to the exper- 
imental results. These calculations are included in Fig. [J] 
(dotted red line). 

There are two ways to interpret the dominant isoscalar 
term along with a nearly zero isovector contribution to 
the B(E2) strengths. Conceptually, this may be per- 
ceived as the two-alpha cluster in these nuclei being the 
same in 10 Bc and 10 C and their respective valence parti- 
cles contributing little. This mirror-symmetric interpre- 
tation conforms well to the predictions of the standard 
shell model. One can alternately consider the possibility 
that mirror symmetry is not preserved; for example, the 
alpha-cluster spacing is modified, but this effect is off- 
set by differing contributions from the valence particles. 
This effect can be investigated through the use of more 
sophisticated models. 
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FIG. 4: (Color online) B(E2; 2+ 0+) transition strengths 
(in e 2 fm 4 ) for 10 C and 10 Be. Open symbols are the results 
of Ref. 01 while solid symbols are the current work and the 
recent measurement of Ref. [J| . 
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The variational Monte Carlo (VMC) and Green's func- 
tion Monte Carlo (GFMC) methods have been very use- 
ful in improving our understanding of light nuclei fl6l fT7| . 
and successful in reproducing the electric quadrupole 
collectivity in 10 Be, without resorting to the use of 
any effective charges. Using realistic two- and three- 
body forces and operators (including explicit chargc- 
symmetry-brcaking terms), this is a good approach for 
exploring 10 C and the symmetry of the A = 10 wave 



functions. More details on the GFMC method of calcu- 
lating transition strengths are given in Ref. [l8j . 

The VMC calculations use trial wave functions con- 
taining non-central, two- and three-body correlation op- 
erators acting on an antisymmetrized one-body wave 
function, $(JA/TT Z ), which determines the quantum 
numbers of the state being computed. The &(JMTT Z ) 
wavefunction is expanded in LS'-basis functions [l7| : 



$(JMTT Z )= P( 2S+1 L[n},JTT z )<£>( 2S+1 L[n},JMTT z ), (2) 

LS[n] 

where the amplitudes /3( 2S+1 L[n], JTT Z ) are found from a diagonalization of the Hamiltonian. For 10 Be, we construct 
states from the three highest spatial symmetries as denoted by the Young diagram [n] (see Ref. [HI). This gives 
three basis functions for the + ground state: 1 S[442], 3 P[4411], and 3 P[433], while the 2 + states have six basis 
functions: X D[442]_, 1 D[442]+, 3 P[4411], 3 P[433], 3 F[4411], and 3 F[433]. Note that there arc two linearly independent 
1 D[442] basis states; the distinction between them is arbitrary and we choose to express them as eigenfunctions of the 
quadrupole operator, with the subscript indicating the sign of the quadrupole moment. The VMC E2 matrix element 
is 

M(E2)= f3( 2S ' +1 L'{n%J = 0,T=l,T z )f3( 2S+1 L{n},J = 2,T=l,T z ) 

L'S'[n'],LS[n] 

x(C$( 2S ' +1 L'[n'},J = 0,T=l,T z )\\E2\\C${ 2S+1 L[n},J = 2,T=l,T z )} , (3) 



where C denotes the two- and three-body correlations and 
there are 3x6 contributions to the sum. The E2 operator 
does not change spatial symmetry, so the only big contri- 
butions are those from !D[442]_ or ^[442]+ to X S[442], 
3 P[4411] or 3 F[4411] to 3 P[4411], and 3 P[433] or 3 F[433] 
to 3 P[433] (the C do not conserve the spatial symmetry 
so there are small non-zero matrix elements for the other 
possibilities). These individual contributions, calcu- 
lated with wave functions for the AV18 two-nucleon and 
Illinois- 7 three-nucleon potentials (AV18+IL7) [H [Hj], 
are shown in Fig.[S]for isospin-symmetric basis states, i.e., 
the parameters in $( 2Sfl £[n], J = 2, M,T=1, T z ) are in- 
dependent of T z . The diagonalization of the two 1 D[442] 
states into the quadrupole basis was done for 10 Be and 
not changed for 10 B and 10 C. As can be seen, these 
isospin-symmetric calculations can give very different T z 
behaviors, depending on the pair of 2S+1 L[n] states being 
used. Calculations with basis states containing different 
variational parameters give very similar results to those 
in the figure; we believe that the trends shown result 
from the different 2S+1 L[n] values of the pairs and thus 
also would be obtained with other realistic Hamiltonians 
and even for the corresponding harmonic-oscillator shell 
model states. This means that the nearly T z -independent 
B(E2) strengths obtained in the shell model calculation 
require a specific combination of states. 

Figure[S]shows VMC and GFMC calculations oiB(E2) 
values for the 2^~ (Q < 0) state assuming isospin symmet- 
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FIG. 5: (Color online) B(E2) transition strengths (in e 2 fm 4 ) 
for VMC basis state pairs. Only transitions that conserve 
spatial symmetry are shown. The legend gives the designation 
of the + states while the curves are labeled with the 2S+1 L 
of the 2 + states. 



ric wave functions. The wave functions were computed 
for 10 Be using the AV18 interaction alone or with the IL2 
or IL7 three-body potentials. As shown in Ref. the 
2^~ state of 10 Be has a negative quadrupole moment and 
a strong E2 decay to the ground state for all the Hamilto- 
nians. (For AV18 alone, the energies of the two 2 + states 
are nearly degenerate, so we choose to identify the Q < 
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FIG. 6: (Color online) VMC and GFMC calculations of 
B(E2) strengths for the 2f (Q < 0) state assuming isospin 
symmetric wave functions. The wave functions were com- 
puted for 10 Be using the indicated Hamiltonians. VMC re- 
sults are shown as open symbols and dashed lines; GFMC 
results are solid symbols and solid lines. 
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FIG. 7: (Color online) VMC and GFMC calculations of 
B(E2) for the lowest (Q < 0) 2+ state with the AV18+IL7 
Hamiltonian. The wave functions were computed separately 
for each nucleus using different VMC correlations and inde- 
pendent GFMC propagation. VMC results are shown as open 
symbols; GFMC results are solid symbols. The dashed lines 
show the corresponding isospin symmetric results from Fig. [6] 



state as the 2f state.) The reduced matrix elements for 
10 C were obtained by interchanging protons and neutrons 
in the 10 Be wave functions and the 10 B reduced matrix 
elements are the average of the 10 Be and 10 C, i.e., the 
wave functions are isospin symmetric with those of 10 Be. 
There is considerable variation in the T z behavior of the 
B(E2) strengths for the different Hamiltonians in the 
VMC calculations. This is presumably due to the differ- 
ent (3( 2S+1 L[n], J=2,T=1,T Z = +1) amplitudes from 
the separate diagonalizations. The GFMC generally pre- 
serves, or even enhances, these different trends which sug- 
gests a strong sensitivity of the isovector B(E2) to the 
three-body force. 



Isospin symmetry of the wave functions is certainly 
only an approximation. Due primarily to the increasing 
Coulomb potential energy going from 10 Be to 10 C, the 
10 C states under consideration are less bound (0.5 MeV 
vs. 3.5 MeV for the 2\ state) and, hence, should be 
more diffuse. This can be studied by performing sepa- 
rate calculations for each nucleus. We have done such 
calculations for the AV18+IL7 Hamiltonian. The one- 
body parts of the VMC wave functions are solutions 
of Woods-Saxon wells plus an average Coulomb poten- 
tial (T3| ; the strength of the Coulomb term is propor- 
tional to the number of p-shell protons. Separate diag- 
onalizations were made for each nucleus, so the /3 are 
also different. The GFMC propagations are still made 
in a good isospin basis, but the isoscalar Coulomb po- 
tential used reflects the total charge of the nucleus (l8j |. 
The results of these calculations are compared with the 
isospin-symmctric AV18+IL7 calculations in Fig. [7J The 
independent calculations for 10 C and 10 B are not very 
different from the isospin-symmetric extrapolations from 
the 10 Be results. Unfortunately, the already too large 
value for the 10 C ~B(E2) is further increased. 

In light of the apparent failure of these GFMC cal- 
culations to reproduce the B(E2; 2^ —> 0^~) transition 
strength in 10 C, it is important to consider possible short- 
comings of the calculation which could cause the discrep- 
ancy. One possibility is that, with the weaker binding of 
the 2f state in 10 C compared to 10 Be, contributions from 
beyond the p-shcll might become important. In fact, the 
VMC trial functions already have a fair admixture of 
sd-shell and higher components due to two-body tensor 
correlations in the C of Eq.[3](see Ref. [H|), and these are 
further enhanced in the GFMC propagation. However, 
to further test this possibility, we constructed alterna- 
tive clustcrized VMC trial functions with explicit sc?-shcll 
components. These wave functions combine a 8 Be(0 + 
or 2 + ) core with two final nuclcons in p- or sd-shcll or- 
bitals (with appropriate Coulomb terms), all iS-coupled 
to give the appropriate total J n . The + (2 + ) states 
have four (seven) p-shell and eight (seven) sd-shell com- 
ponents; separate diagonalizations for the corresponding 
/3s are made for each T z . For a Hamiltonian contain- 
ing AV18 and the Urbana IX (UIX) three-body poten- 
tial [12], the sd-shell /3s contribute only 2.5% of the total 
wave function in the 10 Be 0+ state, but 21% in the 0^" 
state; these numbers increase to 3.4% and 28.9%, respec- 
tively, in 10 C. The 2f and 2^~ states both have 4% or less 
s<i-shell contributions, with only slightly greater amounts 
in 10 C than in 10 Bc. (Interestingly, 96% of the 2f state 
has a 8 Bc(2 + ) core in this construction, clearly indicat- 
ing it is the J=2 member of the K=0 rotational band.) 
Consequently, inclusion of the sd-shell components has 
only a very minor effect on the B(E2) values; in 10 Bc 
using p-shell only components gives B(E2) = 8.6 e 2 fm 4 , 
while adding sd-shell components raises it to 9.2 e 2 fm 4 . 
The corresponding B(E2) values in 10 C are 9.6 and 11.6 
e 2 fm 4 , respectively, showing the same moderate change 
with T z as the AV18+IL2 and AV18+IL7 Hamiltonians. 
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TABLE II: GFMC calculations of the A = 10, T = 1, ground-state energies E gs and excitation energies in MeV for 
several Hamiltonians used in this work, plus NCSM results for the CD-Bonn 2000 potential and experimental values. ( 10 B 
excitations are shown relative to the + ; 1 isobaric analog.) Also shown are ground state charge radii in fm, quadrupole moments 
for the excited states in efm 2 , and the B(E2) transition strengths in e 2 fm 4 . Asterisks denote GFMC results obtained with 
isospin-symmetric wave functions generated from 10 Be (see text). 
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A more likely possibility is simply that the Hamilto- 
nians tested are not adequate for these transitions. The 
first priority in theoretically modeling the nuclear Hamil- 
tonian has been obtaining good energies for the states in 
question - both absolute binding energy of the nucleus 
and excitation energies of the higher states. These en- 
ergies are shown in Table [TT] for the various Hamiltoni- 
ans used in the present GFMC calculations, along with 
charge radii, quadrupole moments, and the B(E2) val- 
ues. The AV18+IL7 Hamiltonian gives a particularly 
good overall reproduction of both absolute binding and 
excitation energies. The charge radius for 10 Be is also in 
excellent agreement with a recent measurement [23j . Ta- 
ble |H] and Fig. [5] indicate that while the various models 
tested give rather similar results for the B(E2) strengths 
in 10 Be, they give a much more widely varying range of 
results for 10 C. The quadrupole moments also have far 
more variation in 10 C. 

In Table ITU we also show results for the CD-Bonn 2000 
(CDB2k) potential u4[ evaluated with the no-core shell 
model (NCSM) ^M^. These B(E2) results are in ex- 
cellent agreement for both the transitions in 10 Be and 
also very good for 10 C, although they note that sepa- 
rating the two 2 + states in 10 C is non-trivial, leading to 
large error bars for the transition and quadrupole mo- 
ment. The ground state energies are significantly under- 



bound, which is not surprising for a two-nucleon force 
alone, but the 2 + states are reasonably well separated. 
The 10 Be charge radius is also too small, probably due 
to the practical limitations of using a finite harmonic- 
oscillator-spacc basis. 

Many other groups are developing new methods for cal- 
culating light nuclei. They range from models based on 
effective field theories (29[ and density functionals [30| . 
to monte carlo shell models [3l| and cluster-based ap- 
proaches (32l . [33j . New, precise experimental measure- 
ments of masses, radii, moments and decays can all 
help to refine these approaches and improve our insight 
into nuclear structure and its evolution with neutron-to- 
proton ratio. 



IV. CONCLUSIONS 

We have performed a precise measurement of the life- 
time of the first (and only) excited bound state in 10 C. 
The new measurement implies a matrix element that is 
considerably smaller than previously reported. It is only 
2% different from its mirror transition in 10 Be. GFMC 
calculations with our best Hamiltonian fail to reproduce 
this near equality. Different basis states have very dif- 
ferent trends of the B(E2) going across this isomultiplet. 
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Thus, if the wave functions are nearly isospin symmet- 
ric, our calculations fail to get the correct mixture of 
states. Inclusion of explicit sd-shell amplitudes into the 
trial wave functions, while slightly different in 10 Be and 
10 C, appears to be a small effect that does not resolve the 
discrepancy. The calculations show only a small amount 
of isospin symmetry breaking. The definite experimental 
signal of such symmetry breaking for this isovector tran- 
sition operator would be an observation that the matrix 
element for the transition between the first T — 1, 2 + 
and + states in 10 B is not the average of those for 10 Be 
and 10 C. Such a measurement is underway. 

Obtaining a better theoretical result would appear to 
require a Hamiltonian that produces a very specific com- 
bination of spatial symmetry components. It may be pos- 
sible to develop improved three-nucleon potentials that 
provide both good energies and transition strengths at 
the same time. Recently, it has been suggested that us- 
ing the weak /3-dccay of tritium in addition to the three- 
nucleon binding energies is a good way to select an opti- 
mal combination of low-energy constants in chiral N3LO 
three-nucleon potentials (28|. In a similar manner, elec- 



tromagnetic transitions like the ones studied here may 
provide a valuable constraint on suitable Hamiltonians, 
i.e. on the 3-body term. 
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